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ABSTRACT 

We investigate the correlation between the mass of a central supermassive black 
hole and the total gravitational mass of the host galaxy (Mtot)- The results are based 
on 43 galaxy-scale strong gravitational lenses from the Sloan Lens ACS (SLACS) Sur- 
vey whose black hole masses were estimated through two scaling relations: the rela- 
tion between black hole mass and Sersic index (Mbh — n) and the relation between 
black hole mass and stellar velocity dispersion (Mbh — c^)- We use the enclosed mass 
within R20O) the radius within which the density profile of the early type galaxy ex- 
ceeds the critical density of the Universe by a factor of 200, determined by gravita- 
tional lens models fitted to HST imaging data, as a tracer of the total gravitational 
mass. The best fit correlation, where Mbh is determined from Mbh ~ relation, is 
log (Mbh) = (8.18 ± 0.11) + (1.55 ± 0.31)(log (M^t) - 13.0) over 2 orders of magnitude 
in Mbh- From a variety of tests, we find that we cannot reliably infer a connection 
between Mbh and Mtot from the Mbh — n relation. The Mbh — Mtot relation provides 
some of the first, direct observational evidence to test the prediction that supermassive 
black hole properties are determined by the halo properties of the host galaxy. 

Subject headings: black hole physics, gravitational lensing, galaxies: evolution, galaxies: 
halos, quasars: general 



1. INTRODUCTION 



Supermassive black ho les (SMBH) are believed to reside in nearly all galaxies (iKormendy &: Richstone 



1995l : lFerrarese et al.ll2005l ) and the masses of these SMBH (Mbh) show correlations wit h host galaxy 



properties, implying that SMBH and ga. 



2001 



2000 



axy formation processes are closely linked (I Adams et al 



Hopkins et al.l i2005a b; ISilk et al 



1998 



CattaneolboOllj ICattan eo et al.|[l99g:lDi Matteo et al.ll2005l.l2003l:lEl-Zant et al.lEo03, : .Haehnelt fc Kauffmann 



Wvithe k, PadmanabhanI l2006l ) . Previous studies 
have shown correlations between Mbh and galaxy's effective stellar velocity disp ersion (a^,), bulge lu- 



minosity (Lbui), Sersic index (n) and stellar mass of the bulge component (Mbui) (jFerrarese h Merritt 
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2000; iGebhardt et alJbood : ICraham et alJl200ll . 120071 : iMarconi et al-lbood ) . Some of the challenges 
faced by current models of SMBH formation and evolution include reproducing and maintaining 
these scaling relations regardless of the ev ents that take place during galaxy evolution driven by 
the process of hierarc hical mass assembly (jCrotonl l2009l : IWvithe et al.l |2002| . l2003l : iMcLurd 120061 : 
Robertson et al.ll2006l ). These scaling relations are not only important tests of the models of SMBH 
formation and evolution but also powerful predictive tools from which we can infer SMBH masses 
of galaxies that are located at higher redshifts. 

In this paper, we examine the evidence for a scaling relation between Mbh and the total mass 



predict a fundamental connection between Mhh and M+ot of the host galaxy ( 


Adams et al. 


2001: 


El-Zant et al. 


2003: 


Haehnelt et al. 


1998: 


Monaco et al. 


2000: 


Silk et al. 


1998) 


. One of the most 



important predictions of the galaxy models, which study the interaction between the dark matter 
haloes of galaxies and baryonic matter settling into the gravitational potential to form the bulge 
and SMBH, is that halo p r operties determine th ose of the bulge component and SMBH (ICattaneo 



2001 



El-Zant et al. 



2003 



Hopkins et al. 



2005 



3H). 



However, observational evidence for such a 



scaling relation has been sparse since measurement of the total mass is non-trivial. 



Results from iFerraresd ()2002l ) examir ied some o f the first indirect observational evidence for 
the existence of a Mbh — Mtot relation. In IFerraresd (120021 ). the correlation between bulge velocity 
dispersion (ctc) and the observed circular velocity (vc,obs)) for a sample of 20 elliptical galaxies and 
16 spiral galaxies, is translated into an equivalent Mbh — Mtot correlation. Although ac can be 
translated into Mbh through Mbh — relation in a straightforwards manner, an estimate of Mtot 
is dependent on the conversion betw een v^^obs and virial velocity, the velocity of the galactic halo 
at the virial radius . IFerraresd (120021 ) uses results from ACDM cosmological simulations to derive 



Following IFerraresd (j2002l ) , several stud ies have examined the Vp^nhs — q"r relation and its 



Mt, 

implications for galaxy fo rmation and evolution (jBaes et al. 
200?! : IPizzella et allbood ). 



2003 



Buvle et al 



2006 



Courteau et al 



The primary goal of this study is to extend the observational evidence to support the Mbh — Mtot 
relation using an attractive, alternative, and direct method to measure the total mass of the host 
galaxy. In recent years, strong gravitational lensing has emerged as a powerful tool to probe the 



mass prof i les, ranging from individual galaxies to clusters of galaxies (iBolton et al 



Dve et al 



2007 



Gavazzi et al 



200S 



Halkola et al 



2008 



Moustakas et al 



2006 



2008a 



20071 ). Due to the fact 



deflection of a photon passing an intervening massive object is independent of the deflector's dy- 
namical state, gravitational lensing does not suffer from difficulties associated with dynamical mass 
measurements of galaxies, where assumptions regarding orbital motions of tracers can lead to 
complications. Thus, gravitational lensing is a more robust method to estimate the total mass, 
including dark matter, around early - type gala.xies and infer the existence of an isothermal mass 
profile {p oc r~^) in various systems ( Dye et al.ll2007l : iKoopmans et al.ll2006l ) 



The mass enclosed within the Einstein radius (Meinst) which is measured from the lens model 
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that produces the best fit to the observed multiple images, is a direct probe of both luminous and 
dark mass in a galaxy. In conjunction with velocity dispersion and surface brightness profiles of 
the lens galaxy, measurement of Meinst from a lens model can be used effectively to constrain the 
luminous and dark matter profiles. In this study, we utilize the Meinst as a tracer of the total 
mass and use the total mass profile of the lensing galaxy to determine the mass contained within 
a redshift-independent circular aperture. 



We also attempt to derive the Mbh — Mtot relation using purely photo metric tracers of M 



and Mtot- Therefore, we examine the possibility of using Mbh — n- relation (jGraham et al.l 12001 



20071 ) as the primary method to estimate black hole masses. Surface brightness profile s of the bulge 
component of most galaxies, in particular E/SO, can be described by the Sersic law (jSersidll968l ) 
as follows: 



i;(r) = Socxp (-k[(r/rc 



1]) 



(1) 



where S(r) is the surface brightness at radius r and Sersic index (n) is the degree of light concen- 
tration. Previous studies in the literature indicate that the quantity n varies monotonically with 
galaxy magnitude. Therefore, the existe nce of a Mhh — n rela tion can be inferred given the depen- 
dence of Mbh on the galaxy magnitude (Marconi et al.l 120031) and the connection between gala xy 



magnitude and Sersic index (jjerien et al 



2000; iGraham k Guzmanll2003l : iFerrarese et al.1 12009 ). 



Graham et al.l (120071 ) discuss the most recent version of the Mbh — n relation and indicate 
that the dependence between the quantities log(Mbh) and log(n) can be best represented by a log- 
quadratic relation. The main motivation for using this relation is that measurement of n requires 
only imaging data, which is easier to acquire than spectroscopic data at a given redshift; therefore, 
Mbh " Mtot could be easily extended to large samples of higher redshift gravitational lenses, which 
may not have spectroscopic measurements. Furthermore, n is a distant-independent quantity, an 
added advantage for estimating black hole masses of high redshift galaxies that may not have secure 
redshift measurements. 

This paper is structured as follows. In §2 we present a brief description of the target selection. 
In §3 we describe the following analysis procedures: §3.1 Measurements of Mtot! §3.2: Deriving the 
connection between Mbh and Mtot using Mbh ~ cr* relation; §3.3: Deriving the connection between 
Mbh and Mtot using Mbh ~ ^ relation. In §4 we present the results of our analysis. In §5 and §6 we 
discuss the implications of our results and give some concluding remarks regarding this project. We 
assume the following cosmological terms for all computations in the paper: J7m = 0.3, = 0.7, 
Hq = 70 hyo kms^^ Mpc^^ and hyo = 1. Unless otherwise noted all scaling relations in this paper 
are defined as linear relationships in log-log space and all logarithms assume a base of 10. 
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THE SAMPLE 



Our sample i s a subset of the ga l axy-sc ale strong gravitational lenses from the Sloan Lens ACS 



(SLACS) Survey (jBolton et al.ll2006l . l2008al ). a Hubble Space Telescope (HST) imaging survey that 
is optimized to detect bright early-type lens galaxies with faint lensed sources. These lenses were 
initially selecte d from spectra of gala.xies of the SPSS Lu minous Red Galaxy (LRG) and MAIN 



galaxy samples (lEisenstein et al.ll200ll : IStrauss et al.ll2002l ). After a best-fit template of the contin- 



uum was subtracted from the spectrum of each galaxy, the residual spectra were scanned for nebular 
emission lines at a redshift higher than that of the target galaxy. Galaxies whose residual spectra 
exhibited at least three atomic transiti ons at a common background redshift were identified as lens 



candidates ([Bolton et al.l l2006l . l2008al ) and imaged using the Wide-Field Channel (WFC) of the 



Advanced Camera for Surveys (ACS) through F435W and F814W filters (Programs 10174, 10587 



subsequent analysis can be found in the series, The Sloan Lens ACS Survey 1-8 i 


Bolton et al. 


2006. 


2008a 


b; 


Gavazzi et al. 


2007. 


2008: 


Koopmans et al. 


2006; 


Treu et al. 


2006. 


2009 


). To date, SLACS 



survey has discovered 131 galaxy-scale strong gravitational lenses with snapshot or deep HST ACS 
imaging, secure spectroscopic redshifts for both lens and source galaxies in ea ch system and stella r 



velocity dispersion, cjsdsS) measured within the 3" SDSS spectroscopic fiber (jBolton et al.ll2008al ). 



From the 131 SLACS lenses, 70 systems are classified as "Grade- A" systems bv lBolton et al 



(|2008al ). where the Grade- A classification implies a definite lens as determined by HST imaging. 
We use 43 "Grade- A" lenses, from programs 10174, 10587 and 10886, as the sample for our analysis. 
We do not use lens systems that do not have a simple lens model due to nearby companion galaxies 
(7 lens systems) or systems where cjsdss value is not available (14 lens syste ms) because the med ian 
signal-to- noise (SNR) was too low within the 3" SDSS spectroscopic fiber (jBolton et al.ll2008al ). 



DATA ANALYSIS 



The most important advantage of using the SLACS galaxy sample to determine the Mbh — Mtot 
relation is that we have an accurate measurement of the total mass profile of each galaxy. However, 
SLACS lenses are located at distances ((ziens) ~ 0.2) where the SMBH mass cannot be measured 
directly through its influence on the surrounding gas or stars. Therefore, we use two scaling relations 



to determine the SMBH masses of the gala-x y sample: M^h — <7* rel ation (IFerrarese Merritt 
Gebhardt etTallboOol : Iciiltekin et, aPbood ) and Mbh - n relation (jCraham et ahlboOll . I2OO7I ). 



200c 



3.1. Deriving Total Mass of the Host Galaxy 



We us e parameters f rom th e strong lens modeling performed by iBolton et al.l (j2008al ) to ob- 
tain Mtot- iBolton et al.l (j2008al ) use a singular isothermal ellipsoid (SIE) lens model to describe 
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foreground lens plane which generates the multiple images of the source. The SIE lens model is 
parameterized by the angular Einstein radius, bsiE; which relates to the physical parameters of the 
mass model as follows: 



bsiE 



47r^ 



Dls 



(2) 



Dos 

where Dls and Dqs are angular diameter distances between the lens and source plane and observer 
and source plane respectively. The surface brightness profile of the background galaxy is described 
by a single or multiple Sersic ellipsoid profile and the model lensed image, for a given SIE lens 
model, is produced by the forward ray tracing method. The model parameters are adjusted until 
the images match the o bservations and the solutions are optimized using a merit function, such as 



(IBolton et al.ll2nn8al ^ 



We use the lens model parameters given bv lBolton et al.l (j2008al ). i.e. bsiEi and integrate the 
best-fit mass profile to obtain the total mass (both luminous and dark) within a projected radii 
of choice. During the preliminary analysis, we use the inverted lensing equation to compute mass 
inside the Einstein radius (Meinst)) as a tracer of Mtotj assuming that the lens and source galaxies 
are aligned along the optical axis. 

However, Meinst is not a measurement of mass within a standard aperture since the projected 
radius characterized by the Einstein radius for a given mass distribution is redshift dependent. We 
therefore use R20O) the radius within which the mean density of the lensing galaxy exceeds the 
critical density of the Universe (pcrit) by a factor of 200, as a "standard aperture". The use of R2nn 
as a radial aperture is also consistent with theoretical studies of the Mbh — Mtot relation (jCroton 
2003 ). We integrate the best-fit mass profile to obtain the mass contained within a projected radius 
of R200 (M200) as a tracer of the total mass as characterized below: 

2 (<7sie)^ R200 



M 



200 



G 



(3) 



Adopting an "aperture-corrected" mass measurement (M200)! instead of Meinst; does not affect 
the overall form of the Mhh — Mt.nt r e lation (the slope and zero-point) discussed in the follow- 



ing sections. Results of iBolton et al.l (|2008bl ) show a lack of correlation between the quantities 
f = o-aporture/o-siE and mass or Reinst/Raperture, where (Taperture is the Stellar Velocity dispersion cor- 
rected to an aperture of Raperturc- These results are consistent with the near isothermal nature of 
the radial profile; thus, using a radial aperture other than Reinst, to derive the total mass, does not 
introduce an additional dependence on the assumed mass model. We also observe a relationship 
of log (Meinst /M200) oc log (Reinst/R20o)> Consistent with the fact that an aperture corrected mass 
measurement does not affect the observed Mbh — Mtot relation. 

To obtain the associated 1-a error of M20O) we use Ga ussian error propaga tion adopting an 
empirical error of 2% on the measured Einstein radii following lBolton et al.l (|2008al ) {{5 log(M2oo)) ~ 
0.01 dex). Table 1 lists the unique SDSS identifiers for SLAGS lenses, redshift of the lens galaxy of 
a system and the derived values of M200 (hereafter denoted as Mtot) for each lensing system. 
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3.2. Estimating SMBH Masses Using Stellar Velocity Dispersion 



We first use the Mbh — o"* relation to estimate Mbh of the early type lens sample. The velocity 
dispersion measurements are derived from spectra from the SDSS 3" diameter fiber. The SDSS 
velocity dispersion measu rements are corrected t o a standard aperture, of radius equal to re/8, using 
the power-law relation bvlJorgensen et al.l (119951 ). similar to SMBH and Fundamental Plane studies 
(jFerrarese et al.l l2005l : iTreu et al.l 120061 ). The mean difference between SDSS velocity dispersion 
measurements and the aperture corrected measurements (hereafter denoted as o"^) are ~ 4%; 
therefore, the aperture correction does not significantly influence the overall form of the Mhh —_M^t 
relation discussed in the following sections. We use Mbh — relation by iGiiltekin et al.l (j2009l ). 
obtained from measurements of Mbh from dynamically detected central black holes: 



log(Mbh/M0) = (8.12 ± 0.08) + (4.24 ± 0.41) log(o-^/200 km s 



(4) 



with an intrinsic scatter of Eq = 0.44 dex. 

We assume that errors on cj^ and best-fit coefficients of the Mbh ~ cr* relation are uncorrelated 
and use Gaussian error propagation to determine the 1-a error on the quantity log (Mbh)- The values 
of (Tv, and its associated error for each lensing system and SMBH mass estimates obtained from 
Mbh ~ o"* relation are listed in Table 1. 



3.3. Estimating SMBH Masses Using the Sersic Index 



We also attempt to u se the Mbh — relation to estimate Mbh of the SLAGS lens sample 
(jGraham et al.ll200ll . 120071 ). We perform a careful and detailed decomposition of HST AGS images 
to obtain the Sersic indices of the SLAGS lens galaxies. Although our sample of lens galaxies were 
observed in at least two ACS filters, mainly F814W and F435W, we exclusively use F814W (Broad 
I-band) data due to high SNR of the images. Furthermore, redder magnitudes are a better tracer of 
m ass in comparis o n to B -band magnitudes. To keep our decomposition method consistent with that 
of lGraham et al.l (|2007l ). we use a two-component model to describe the surface brightness profiles 
of SLAGS lens galaxies. In addition to the Sersic profile, which describes the surface brightness of 
the bulge component in each galaxy, we include a second disk component characterized by a simple 
exponential profile: 



S(r) = Soexp(-r/rd) 
where E(r) is the surface brightness at radius r. 



(5) 



To obtain the best possible accuracy of the data analysis process, we take particular care 
during all intermediate steps leading to the bulge-disk composition. Therefore, we obtain the raw 
images and perform careful cosmic ray removal, distortion correction, manual mask production and 
determination of the best possible AGS point-spread function (PSF). We obtain the AGS images 
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from the HST archive at The Canadian Astronomy Data Centre (CADC), where the images are 
processed by CALACS caHbration software. Imaging data from programs 10174 and 10587 are 
420 s single-exposure snapshot data; therefore, we perform an addi tional cosmic ray identification 
step using the L. A. Cosmic software (LACOS) ( van Dokkum 200ll ). which distinguishes between 



undersampled sources and cosmic rays. We then use the MULTIDRIZZLE reduction package to 
obtain distortion-free ACS images, where the distortion solution is applied to correct all pixels to 
equal areas. 



We use the photometry package SExtractor (jBertin et al.lll996l ) to produce a catalog of galaxies 



in each ACS image. Furthermore, we use the segmentation images generated by SExtractor, which 
deblends each source in the field, to produce manual masks of the lensed features around the early 
type galaxy before lens modeling. We perform the two-dimensional decomposition of e ach early- type 



lens galaxy using Galaxy IMage 2D (GIM2D), which is publicly available to users (ISimard et al 



20021 ) . GIM2D uses the Metropolis Algorithm to derive the best-fit parameter values and confidence 
intervals, for a surface brightness model of choice, through Monte-Carlo sampling of the likelihood 
function. Since our surface brightness model consist of a bulge and a disk component, we explore 
the full range of bulge fraction (B/T = 0.0 - 1.0) such that the Metropolis Algorithm converges 
to an accurate quantitative morphology of galaxies classified as E-types. On average, we find that 
most SLACS lenses have a significant bulge component (( B/T ) 0.73) and a bulge plus disk light 
profile produces a better fit to the lens galaxy, significantly minimizing the residuals when the 
galaxy model is subtracted from the input galaxy image. 

For GIM2D lens models, we take particular care to define the ACS Point Spread Function 
(PSF), which is extremely position- and time- dependent. We investigate a variety of PSF models 
to determine the most suitable PSF for convolution with GIM2D galaxy models. The PSF models 
are as follows. 



1. A star extracted from the field. 



A PSF p roduced at the location of the galaxy, in the non-drizzled image, using Tiny Tim 
software (lKristill993l ). We insert the Tiny Tim generated PSF into an empty ACS-WFC field, 
at the location of the galaxy, and multidrizzle the resulting image to obtain a distortion-free 
PSF model. 



3. A PSF produced at the location of the galaxy u sing the principle component analysis (PCA) 
models in the ACS PSF library (jjee et ahlboOTl ). 



We examined the residual images produced by subtracting the lens surface brightness model, each 
convolved with the PSF models described above, to determine the best-fit PSF for GIM2D lens 
modeling. The PSF model produced from the ACS library gives the best results, minimizing the 
core residuals for a fixed surface brightness model. Therefore, we extract the Sersic indices from 
the best-fit bulge plus disk decomposition, using the ACS PSF library model for convolution. An 
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example of the difference between residual images produced by various PSF models are shown in 
Figure [H 

We also examine the reliability of the bulge-disk decomposition by comparing our results for 
the SLACS lens sample relative to the scaling relations for spheroidal components of the SDSS 
galaxies. Since SLACS lenses were initially derived from SDSS Luminous Red Galaxy (LRG) 
and MAIN samples, structural parameters measured from ACS imaging should be in agreement 
with those of SDSS galaxies. Figure [2] shows the comparison of effective bulge radius (Re), Lband 
magnitude of the bulge component (Mi^buige) and aperture-corrected bulge velocity dispersion (cTap) 
of the SLACS lenses to the SDSS bulge parameters. We construct the scaling relations for SDSS 
galaxies from the bulge-disk decomposition of 77523 galaxies and velocity dispersions from the 
SDSS database. The structural parameters derived from ACS imaging are in good agreement with 
those of SDSS galaxies, which confirms the accuracy of our bulge-disk decomposition (in addition 
to the minimal residuals seen in the galaxy images after subtracting the galaxy models). The 
outlier in Figure [2] is SDSS J0959 + 0410, a disk galaxy with a small bulge fraction (B/T = 0.14). 
Previous studies, which examine the positions of bulges with disk-like features (also referred to as 
pseudobulges) in structural paramet er space find that pse udobulges with low B/T deviate from the 
scaling relations of classical bulges (jFisher Drorvll2008l ). 



Using the best-fit parameter for n from GIM2D model s, we estimate the SMBH masses using 
the log-quadratic Mbh — n relation by (jCraham et al.l 120071 ): 



log(Mbh/M0) = (7.98 ± 0.09) + (3.70 ± 0.46) log(n/3) - (3.10 ± 0.84)[log(n/3)]^ 
with an intrinsic scatter of Eq = 0.18 dex. 



(6) 



Similar to lGraham et al.l (120071 ) . we assume a measurement error of 20% on the values of n to 
obtain the 1-cr error on log(Mbh)- Table 1 lists the Sersic indices and bulge fractions of the SLACS 
lenses from GIM2D modeling, estimates of Mbh for each lens system. 



RESULTS 



4.1. Mbh,a. - Mtot Relation 



In the following section, we combine the correlation between the primary observable quantities 
((T* — Mtot) with Mbh — relation to derive the Mbh — Mtot relation. A tight correlation between 
the quantities log((7*/200 kms~^) and log(Mtot) is apparent from Figure El From Spearman's rank 
test, we obtain a correlation coefficient (rg) of 0.84 (degrees of freedom = 41), which indicates a 
strong positive correlation between the quantities log((Tv,/200kms~^) and log(Mtot), with a 99.99% 
confidence level that the correlation has not occurred by chance. 



To quantify this correlation, we use the x^-fitting routine by Weiner et al. ( 2006), which im 



plements a generalized form of the least-squares fitting routine by lPress et al.l (119921 ). 



Weiner et al 



- 9 - 



(|2006l ) routine accounts for intrinsic scatter (eg) beyond the observational errors, for a relation 
of interest, by adding in quadrature to the error in the dependent variable. Initial fits to the 
observed relation between log((T^/200 km s~^) and log(Mtot)) incorporating observational errors in 
both variables, gives a large reduced value (xred^ ~ 3.0) indicative of intrinsic scatter in the 
relation. Therefore, we perform the fits by incorporating observational errors in both variables and 
intrinsic scatter in the y-variable. The value of eg is determined by requiring that Xred^ is unity. 
The results of the fitting routine gives the following best fit correlation: 

log (cj,,/200kms"^) = (0.014 ± 0.013) + (0.365 ± 0.038) [log (Mtot/Mo) - 13.0] (7) 

with an intrinsic scatter of 0.037 dex in log((Tj,/200 km s~^). Combining equations S] and [3 we 
derive the following Mbh — Mtot relation: 

log (Mbh/M©) = (8.18 ± 0.11) + (1.55 ± 0.31) [log (Mtot/M©) - 13.0] (8) 

Figure H] shows o"^ transformed into SMBH masses using equation HI We assume that measure- 
ment errors of the velocity dispersions and best-fit coefficients of the Mbh ~ relation are uncorre- 
lated and use Gaussian error propagation to determine the 1-a error on the quantity log (Mbh)- We 
also perform a direct fitting step to Mtot and the secondary observable quantity, Mbh, to confirm the 
validity of the derived Mbh - Mtot relation. We find log(Mbh/M0) = a + /3 (log (Mtot /M©) - 13.0) 
with (a,/3) = (8.17ib0.13, 1.57±0.39) and Xrod^ ~ 0.2, which is in agreement with the result shown 
in equation El 



4.2. Mbhn - Mtot Relation 



In this section, we discuss the n — Mtot relation of the SLACS lenses , shown in Figure [5l 
and Figure [6] which shows n converted to SMBH masses using equation [6] ([Graham et al.l 120071 ) . 
Similar to §4.1, we use Gaussian error propagation (assuming a measurement error of 20% on n) 
with intrinsic scatter of the Mbh — Mtot relation (eo = 0.18 dex) add e d in q uadrature, to determine 
the 1-0" error of the quantity log(Mbh)- Results of iGraham et al.l (|2007l ) indicate that Mbh — "f^ 
relation is comparable to the Mbh — cr* relation; therefore, we expect both scaling relations to yield 
Mbh — Mtot relations with similar level of scatter. If the expected equivalence between Mbh — n and 
Mbh — fit relations hold, we can then extend the Mbh — Mtot relation to gravitational lens samples 
which span a higher dynamical range in the total mass at various redshifts. Unfortunately, however, 
the trends observed in Figures [5] and [6] are very different from the strong positive correlations shown 
in Figures [3] and [H 

From a Spearman's rank test, we obtain a weak negative n — Mtot correlation (rg ~ —0.42 for 
degrees of freedom = 41) for the trend observed in Figure [5] with only a 95% confidence level that 
the correlation has not occurred by chance. A Spearman's rank test performed on the Mbh — Mtot 
relation shows similar statistical properties to that of the n — Mtot relation. Before concluding if the 
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correlation observed in Figure U] is a better representation of a Mbh — Mtot relation, we investigate 
various possibilities that may cause the inconsistency observed between Figures H] and [6l 



First, we attempted to reproduce the Mbh — n relation, obtained by iGraham et al.l (|2007l ). 
for the SLAGS lens sample. In Figure El we plot Mbh.o-* values for the SLAGS lenses, versus the 
co rresponding Sersic ind ices from the bulge-disk decomposition. The log-quadratic Mbh — n relation 
by iGraham et al.l (|2007l ) is overlaid on the SLAGS data for comparison. Figure [7] clearly indicates 
that the log-quadratic Mbh — n relation cannot be reproduced with the Mbh.o-* and n values for 
the SLAGS lens sample. A weak but different correlation between Mbh.o-* and n is apparent from 
Figure [71 We perform a Spearman's rank test and derive a weak, negative correlation, rg ~ —0.44 
at the 5% significance level, between these quantities. 

The correlation observed i n Figure |7| is in contr ast to the log-linear or log-quadratic Mbh — n 
relation discussed in literature (jGraham et al.l 120071 ). Previous studies of morphology and scaling 
relations of bulge dominated galaxies claim that n scales monot onically with the galaxy bulge 
magnitude (equivalently bulge luminosity) (jFerrarese et al.l 120061 ). If this correlation holds, we 
expect n to correlate posit i vely with Mbh due th e dependence between Mbh and bulge luminosity 
(Lbui) (jMarconi et al.ll2003l : iFerrarese et al.ll2006l ). 



Since Figure [71 contradicts the existence of a log-quadratic Mbh — n relation, we also test 
the correlation between n and I-band bulge magnitude (Mbuige,i) for the SLAGS lens sample. The 
results of GIM2D bulge-disk decomposition of SLAGS galaxies are shown in FigureO A Spearman's 
rank test to determine degree of correlation between Mbuige,i and n shows that there is only a very 
weak correlation between the two quantities (rg ~ —0.20); furthermore, the rank correlation we 
observe between Mbuige,i and n is significant at a level larger than 5%. 

Although we do not observe a direct correlation between galaxy bulge magnitude and the 
Sersic index, our results are limited by the narrow range of magnitudes (—24 < Mbn igej < —20) 
i n the SLAGS lens sample. Furthermore, our results are consistent with those of iBolton et t 
( 2008bl ). who also find that n for all SLAGS lenses is uncorrelated with the measurements of 
lensing mass, dynamical mass, luminosity and velocity dispersion of the sample. Given the trend 
observe d in Figure [HI we cannot co nfidently expect Mbh to correlate with n as claimed by previous 
studies (jGraham et al.ll200ll . l2007l ). The possible origin of Mbh ~ relation has not been examined 
theoretically and is beyond the scope of this paper. For the remainder of this study, we use the 
Mbh — cr* relation to determine the black hole masses. 



5. TOWARDS A Mbh - Mtot RELATION 
The correlation ob served between Mbh and Mtnt i in Figure [H is in remarkable agreement with 



theoretical predictions (iGroton et al 



and local observations (Ferrarese 



Groton 



200 



20061 

"Jof the Mbh 



200S 



Silk et al.lll998l : IWvithe et al 



2002 



20031 ) 



Mtot relation. In the following section, we 
compare our results to several black hole formation scenarios posed by various theoretical studies 
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and examine the implications of our findings. Equation [HI wliicli quantifies the trend observed in 
Figure [U shows that Mbh scales non-linearly with Mtot and that efficiency of black hole formation 
increases with total mas s. Theoretical models that reprod u ce the observed lumi r iosity function 



of hiffh-redshift quasars (lAdams et al.l l2003l : ICattaned I2OOII : iHopkins et al.l l2005al : ISpringel et al, 



2005bl : IVolonteri et al.l l2003l ) predict that Mbh scales as a power law of the circular velocity of the 
galactic halo (denoted as Vc^haio and also referred to as virial velocity) in which the black hole 
resides: 



Mbh OC Vchalo'' 



(9) 



The Mbh ~ Vc^haio relation shown above can be converted into an equivalent Mbh ~ Mtot relation 
by considering the dependence between Vc^haio and halo mass (equivalent to Mtot in this study), 
Vc,haio OC Mhaio^''^- The resulting correlation between Mbh and Mhaio is as follows: 



Mbh CX Mhalo^/' 



(10) 



The circular velocity of a given halo mass is redshift-dependent; therefore, an imp ortant aspect 
of this analytical predict ion is the evolution of the Mbh — Mhaio relation with time (|Crotonll2009l : 



Wvithe et al. 



2002 



2003). 



The slope of the Mbh — Mhaio relation, 7/8, is a valuable indicator of various formation 
scenarios which result in observed black hole populations. In the process of hierarchical mass 
assembly, formation of SMBH is driven by mergers of galaxy haloes. A linear relation between 
Mbh and Mhaio where 7 = 8, results from a formation scenario where the black holes residing in 
the merging haloes coalesce without additional gas accretion. A slope of 7 > 3 is characteristic 
of a merger where the growth of the resulting black hole is dominated by an accretion process 



( Cattaneo 


2001 




Cattaneo et al. 


1999 


• 


Di Matteo et al. 


2005 




2003 


: 


Haehnelt & Kauffmann 


2000: 


Haehnelt et al.l 


19981: 


Wvithe et al. 2002. 


2OO3II. Within this formation scenario, the fraction of 



baryons accreted on to the central regions that feed the black hole is significantly larger for more 
massive haloes, due to a deeper potential well; therefore, massive haloes host larger SMBH. 



In Figure [9l w e plot two analytical predictions for the Mbh — Mtot relation (jCrotonl l2009l : 
Wvithe et al.l[2003 ). which correspond to two formation scenarios, for three different epochs (z = 0.0, 1.0, 5.0). 
Both types of evolutionary tracks indicate that Mbh increases with the total mass of the host galaxy 
at any given redshift and that a galaxy of given total mass hosts a more massive black hole at a 
higher redshift relative to a lower redshift. Within the context of hierarchical mass assembly in 
ACDM cosmology, the decrease in the growth of Mbh relative to Mtot at lower redshifts can be 
caused by processes such as a decrease in merger rates and the gas fraction that is available to fuel 
the central SMBH. 



Another important physical pro cess involved in the evolution of Mbh — Mtnt relation is the 
feed-back regulated growth of SMBH (jpi Matteo et ahlbood : bopkins et alJboOSbl : Isilk et al.l[l99a : 
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Springel et alj l2005al ). A black hole shines at a fraction, rj, of its Eddington luminosity (L^dd) 
following a merger and returns a fraction of this energy into the surrounding galactic gas. A black 
hole shining at its limiting LEdd can unbind the surrounding galactic gas if the energy liberated 
from the black hole is sufficient to overcome the gravitational binding energy of the gas. As the 
mass of a black hole increases through merger driven processes discussed above, the energy output 
can approach the limit where it is sufficient to unbind the entire galactic gas. The unbound galactic 
gas that escapes into the halo is heated beyond the virial temperature and cannot cool during the 
dynamical time of the quasar; therefore, this mechanism can eventually terminate the accretion 
process that feeds the central black hole. 

The growth of a SMBH via a merger-driven, feed-back regulated mechanism implies a relation 
of Mbh oc Vc^haio^ (the req uired rate of energy dep osition to unbind a self-gravitating system is 
proportional to Vc,haio^/G) (iWyithe et al.ll2002l . l2003l ) . leading to a slope of ^ 1.67 in th e Mbh — M^nt 



relation. Evolutionary tracks shown in solid blue lines in Figure [9l (I Wyithe et al.ll2003l ) are examples 
of the formation scenario described above. The evolution ary tracks with a shallower slope, similar 
to the dashed lines shown in Figure [9] (jCrotonI l2009l ) , may be indicative of a modified feed- 
back regulated growth mechanism. If the galactic gas heated by the energy output of the black 
hole cools before the black hole reaches its critical Eddington limit, additional energy (up to a 
factor of c/vc^haio) would be required to unbind the cool gas component. This formation scenario 
leads to dependence of Mhh oc Vp^hain ^, resulting in a slope of « 1.33 in t he Mhh — y^tnt relation. 



It is worth noting that ICrotonI (|2009l ) assumes the Mbh ~ o"* relation by iTremaine et al.l (|2002l ) , 
Mhh oc a^''^^, which yields a relation of Mbh oc Mtot^ '^^- Therefore, we recompute the re sults of 
CrotonI (j2009l ) using the Mbh — o"jr relation used throughout this study (jGiiltekin et al.ll2009l ). which 



gives a dependence of Mbh oc Mtot^'^^- 

Within the context of formation and evolution of SMBH, it is extremely important to accu- 
rately determine the parameters that quantify the slope and the evolution of Mbh ~ Mtot relation. 
These parameters provide significant insight into the dominant formation scenarios that lead to the 
observed black hole populations at various redshifts. Due to the narrow range of redshifts in the 
SLACS lens sample, the Mbh — Mtot relation we derive can be considered as an evolutionary track 
for the mean redshift of the SLACS sample (( z ) ~ 0.2), as shown in Figure [9] (solid black line). 
The Mbh — Mtot relation we derive from the SLACS lens sample provides a unique opportunity 
to compare the theoretical tracks with observational evidence of the Mbh — Mtot relation. This is 
shown in Figure \W\ where we compare the observational Mbh ~ Mtot relation to the theoretical 
predictions of the Mbh ~ Mtot relation at z = 0.2. The dashed region indicates the upper and lower 
1-0" bounds of the observational form of the Mbh ~ Mtot relation found in this study. Inspection of 
Figure [10] shows that our results are in excellent agreement with analytical predictions of feed-back 
regulated growth. However, from our results, it is difficult to distinguish the importance of gas 
cooling in the black hole formation process. 



The Mbh — Mtot relation found in this study is strongly suggestive that halo properties deter- 
mine those of the galaxy and its black hole. This link is also observed in the Fundamental Plane 
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(FP), the two-dimensional projection of the three-dimensional space defined by the quantities of 
surf ace brightness (Ip), effectiv e radius (Re) and central velocity dispersion (o"e), of early-type galax- 
ies (jPiorgovski &: David 119871 ). Results of iBolton et al.l (|2008bl ) show that a sample of 53 SLACS 
lenses define a FP that is consistent with the general population of early- type galaxies from SDSS: 



Ro 



-0.77 



(11) 



The ratio f = fJe/fTsiE for the SLACS lenses is f = 1.1019 ± 0.008 (jBolton et al.ll2008bl : iTreu et al 



20061 ). indicative of a universal isothermal mass profile (also known as the "bulge-halo conspir- 
acy"). Therefore, we can replace the quantities Re and cie in the FP with R200 (given the trend 
observed between the quantities log(Meinst/M2oo) and log(Reinst/R20o)) and itsie- Combining this 
with equations [3] and [U we derive: 



Mbh CX Mtot'-'° le 



(12) 



Ignoring the dependence of equation [12] on le, since it only varies weakly relative to the other 
variables, we can extract the connection between Mbh and Mtot from th e FP of SLACS lens es. 
The slope of 1.30 of the Mbh - Mtot relation derived from the SLACS FP (|Bolton et al.ll2008bl ) is 
in agreement with our results within the 1-a bounds. The manifestation of a Mbh ~ Mtot relation 
within the FP of the SLACS lenses further strengthens the existence of the observed Mbh ~ Mtot 
relation. We require direct measurements of both Mbh and Mtot to fully compare the manifesta- 
tion of Mbh — Mtot relation within the FP and its implications for regularity of early-type galaxy 
formation scenarios. 

An additional advantage of the results from this study lies within the method used to derive 
Mtot • The traditional method of estimating total gravitational mass of a galaxy, used in black hole 
studies, is to convert the observed circular velocity (from galaxy rotation curves) into the velocity 
of the galactic halo (also referred to as the virial velocity). The inferred total mass of a galaxy 
differs depending on the method used to relate observed circular velocity to the virial velocity. The 
effect on the Mbh — Mtot relation from varying assumptions o f the comi e ction between observed 
circular velocity (vc,obs) and virial velocity (vvir), discussed in iFerraresd (120021). is a l so sho wn in 
Figure m The dashed line sh ows the best-fit Mbh — Mtot relation obtained by IFerraresd (j2002l ) using 
cosmological prescriptions of lBullock et al.l (j200ll ) to relate 

Vc.obs and Vvir (Mbh oc Mtot^'^^) and the 
dot-dashed line show s the resu lt ing r elation where Vc^obs = 1.8vvir (Mbh Mtot^'*^)- From the 
relations obtained by IFerraresd (j2002l ). it is evident that the slope of the Mbh — Mtot is affected 
by the method assumed to relate Vc.obs and Vvir- Gravitational lens modeling is independent of 
such dynamical assumptions and provides an elegant alternative method to determine the total 
mass. Therefore, the Mbh ~ Mtot relation characterized by this independent method is not only a 
complementary comparison to existing observational evidence but also provides valuable insight to 
determine the dominant physical processes of SMBH growth. 
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6. CONCLUSIONS 

We use a sample of 43 early-type galaxies, which exhibit galaxy-scale strong gravitational 
lensing, to derive the scaling relation between black hole mass, Mbh, and the total mass, Mtotj of 
the host galaxy. In this study, we use gravitational lens modeling to directly measure the total 
mass and the mass profile of the galaxy rather than converting the observed circular velocity into a 
total mass, the traditional method that is used in black hole studies. We use two alternative scaling 
relations, Mbh ~ cr* and Mbh — ^ to estimate the black hole masses of the lens sample. We obtain 
a tight correlation between and Mtot in the log-log space. In conjunction with the Mbh — cr* 
relation, we derive the observational form of the Mbh — Mtot relation that is consistent with no 
intrinsic scatter. We do not find a significant correlation between n and Mtot- From a variety of 
tests, we find that we cannot confidently estimate black hole masses with the Mbh ~^ relation. The 
scaling relation we observe between Mbh and Mtot is non-linear and is in agreement with theoretical 
predictions of the growth of black holes and observational studies of the local Mbh — Mtot relation. 
The observed Mbh — Mtot relation is also consistent with the studies of the Fundamental Plane of 



SLAGS lenses (iBolton et al.ll2008bl ). which is suggestive of a unified scenario where the properties 
of the host halo determine those of the resulting galaxy and black hole formed through hierarchical 
merging. The observed non-linear correlation between Mbh and Mtot indicates that massive halos 
are more efficient in forming black holes and the slope of the Mbh — Mtot relation is suggestive of 
a merger-driven, feed-back regulated process for the growth of black holes. 

We wish to thank Laura Ferrarese, Chien Peng, Christopher Pritchet, Jon Willis, and the 
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Canada with the support of the Canadian Space Agency. The authors gratefully acknowledge the 
support from the National Research Council of Canada and NSERC through Discovery grants. 
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Fig. 1. — Comparison of varying PSFs on GIM2D galaxy models. Panel 1 shows the HST ACS 
image of an example SLACS lens. The contrast of panel 1 is exaggerated to show the lensing 
features around the early-type galaxy. Panels 2, 3, and 4 show the residual images after GIM2D 
surface brightness model (Sersic bulge + exponential disk) convolved with a multidrizzled Tiny- 
Tim PSF (panel 2), natural PSF (panel 3) and ACS-hbrary generated PSF (panel 4) is subtracted 
from the observed image. Residual image where the galaxy model is convolved with a Tiny-Tim 
PSF shows artifacts such as boxy-core feature, due to the finite size of the Tiny-Tim PSF, evident 
in panel 2. Statistics of the pixels in panel 2 show that the mean value of the core residuals is 
comparable to the background level; therefore, the artifact introduced by the finite Tiny-Tim PSF 
does not significantly affect the quality of the overall fit. In general, we find that both natural 
and ACS-library generated PSFs produce minimal core residuals. However, most SLACS fields are 
relatively devoid of stars, suitable for convolution; therefore, we use ACS-library generated PSF for 
GIM2D galaxy modeling. 
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Fig. 2. — Comparison of SLAGS lenses to the scaling relations for bulge components of SDSS 
galaxies. The red circles indicate the best-fit values for the following structural parameters of the 
SLAGS lenses: effective bulge radius (Re), Lband magnitude of the bulge (Mj^buige) and aperture 
corrected bulge velocity dispersion (a). The scaling relations between bulge structural parameters 
for SDSS early type galaxies are constructed from a sample 77523 galaxies and are shown in gray- 
scale. The gray-scale is scaled to 20% of the peak value of the central distribution to show that the 
structural parameters of the SLAGS lenses lie well within the distribution that encompasses the 
SDSS galaxies. The outlier in this figure is SDSS J0959 + 0410. 



- 20 - 




0.8 - 



5x10 



12 



10 



13 



5x10 



13 



10 



14 



0> 



Fig. 3. — The correlation between cr^ and Mtot) for the sample of early type SLACS lenses. Velocity 
dispersion values from SPSS pipel i ne ar e corrected to an apertm'e of radius equal to re/8 using the 
empirical relation I Jorgensen et al.l (jl995l ). The error bars correspond to the 1-a error of the aperture 
corrected velocity dispersion, from Gaussian error propagation, taking the measurement errors of 
SDSS velocity dispersions in to account. T o tal ma sses of the SLACS lenses are derived from strong 
lens modeling parameters of lBolton et al.l ()2008al ). The l-a errors of Mtoti which are smaller than 
the data points, are incorporated into the fitting routines but not shown in this plot. The solid line 
correspond to the best-fit — Mtot relation for all SLACS lenses in the logarithmic space. 
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Fig. 4. — Same as Figure [3] with cr^ converted into Mbh using Mbh — Mtot relation (jGiiltekin et al 



20091 ). The error bars correspond to 1-a error on the quantity log(Mbh)- The soHd hne corresponds 
to the Mbh — Mtot relation for all SLAGS lenses derived using the cj^ — Mtot and Mbh — M t nt rela - 
tions. The dashed line represents the resulting Mbh — M tpt relation obtained bv iFerraresd (|2002l ). 



where Mtot is computed from cosmological simulations (IBullock et al. 



20011) relating the the ob- 
The dotted line shows 



served circular velocity (fc,obs) to the virial velocity of the host halo {v^ 
the resulting Mbh — Mtot relation if Mtot is computed using fc,obs = 1-8 Vyir (IFerraresd l2002l ) 
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Fig. 5. — The correlation between n and Mtot; for the sample of early type SLACS lenses. The 
Sersic indices are derived from the best-fit bulge+disk decomposition of the SLACS lenses. The 
error bars correspond to the measurement errors on the qu antity n. Total mass es of the SLACS 
lenses are derived from strong lens modeling parameters of iBolton et al.l (|2008al l. The 1-a errors 
of Mtot 5 which are smaller than the data points, are incorporated into the fitting routines but not 
shown in this plot. 
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Fig. 6. — Same as Figure [5] with n converted into Mbh using Mbh — n relation (jGraham et alJl2007l ). 
The error bars correspond to 1-a error on the quantity log(Mbh)-The l-a error of Mtot is smaher 
than the data points, and is not shown in the plot. 
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Fig. 7. — Comparison of log-quadratic Mbh — n relation to SLACS data. Mbh values are derived 
from the Mbh ~ relation and Sersic indices are best-fit results from the GIM2D bulge-disk de- 
composition. Error bar on the legend represent the typical 1-a uncertainty of log (Mbh), derived 
from Gaussian error propagation using measurement errors from the SDSS velocity dispersions 
{{6 log (Mb h)) ~ 0-42 dex) , and a measurement error of 20% on the quantity n. The log-quadratic 
relation by iGraham et al.l (120071 ) is shown by the solid black line. 
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Fig. 8. — The correlation between Mbuige,i versus n for SLACS lens sample. The values of n are the 
best fit bulge Sersic indices from the bulge+disk decomposition using GIM2D. The values of Mbuigc,l 
are the best-fit, extinction-corrected, I-band bulge magnitudes from the bulge-l-disk decomposition. 
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Fig. 9. — The evolution of the Mbh — Mtot relation from vario us analytical p redict ions for 
z = 0.0,1.0,5.0. The pale blue, solid lines are pred ictions f rom IWvithe et al.l (|2003l ). where 
Mbh oc Mtot^'^^- T he dashed lines are predictions from Croton ( 20091 ). where Mv,h oc Mtot^'^^, for 
the same epochs a s Wvithe et al. (j2003l ) (epochs increasing in the same order as shown in the labels 
for IWvithe eTahl (120031 ^). The thick, solid black line is a comparison of the observed Mbh — Mtot 
relation from the SLACS lens sample (Equation [8]) , at ( z ) 0.2. 
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Fig. 10. — Comparison of the observational Mbh — Mtot relation to theoretical predictions of the 
Mbh — Mtot relation at z = 0.2, the mean redshift of the SLACS lens sample. The solid black line 
correspond to the best-fit result from this study and the shaded region shows the 1-a upper and 
lower limits o f the Mhh — M tnt. re lation. The so l id red and blue lines correspond to the theoretical 
predictions of ICrotonI ((20091) and IWvithe et al.l (|2003l ) of the M^h — Mtot relation respectively. 
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Note. — Col. (1): The unique SDSS spectrum identifier for the lens system. Col (2): Redshift of the early-type lens galaxy 
teolton et al.|[2008al) Col. (3): Bulge fraction of the early-type galaxy from the best-fit bulge-disk decomposition. Col (4): 
Sersic index from the best-fit bulge-disk galaxy decomposition. Col. (5): black hole mass estimated from the Mbh — " relation. 
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The 1-(T error on the quantity log (Mbh.n) is a function of n and its measurement errors and can be computed as follows 

jGraham et al.ll2007l ): 51og(Mbh) = ^[log(n/3)]^ + ''"Bt"/^)]'^ ^^^^2 + t^'™-^'^" 'g^;;/^)]^ (^i lAOi + o.is2 Col. (6): Velocity 

dispersion obtained from the SDSS database corrected to a standard aperture of re/8. Col. (7): black hole mass estimated 
from the M^h — a relation. The l-cr error on the quantity log (Mbh.o-*) is a function of ct* and its measurement errors and can 
computed as follows: (5(log{Mbh)) = + (log{o-*/200 km s-i))2 (0.412) _,_ (5 log((T»/200 km s-i))2 (4.242) _,_ q 442 ^ol. 

(8): Total mass of the lensing ga laxy within the pro je cted radius of R20O1 where mass profile of the lensing galaxy is derived 
from the lens modeling results of lBolton et al ] l|2008ah . The typical l-cr error on the quantity log (Mtot) is ~ 0.01 dex. 



